SCIENTIFIC 

REPORTS 




SUBJECT AREAS: 
SYSTEMS BIOLOGY 
BIOPHYSICS 
CELL SIGNALLING 
CELL GROWTH 

Received 
30 January 2013 

Accepted 
8 March 2013 

Published 
26 March 2013 



Correspondence and 
requests for materials 
should be addressed to 
K.A. (k-aoki@lif.kyoto- 
u. ac.jp) 



A Quantitative Model of ERK MAP Kinase 
Phosphorylation in Crowded Media 

Kazuhiro Aoki' 2 , Koichi Takahashi 3,4 , Kazunari Kaizu 3 & Michiyuki Matsuda 1,5 

1 Laboratory of Bioimaging and Cell Signaling, Graduate School of Biostudies, Kyoto University, Sakyo-ku, Kyoto 606-8501 Japan, 
2 PRESTO, Japan Science and Technology Agency (JST], 4-1 -8 Honcho Kawaguchi, Saitama 332-001 2, Japan, laboratory for 
Biochemical Simulation, RIKEN Quantitative Biology Center, 6-2-3 Furuedai, Suita, Osaka, 565-0874, Japan, institute for 
Advanced Biosciences, Keio University, 5322 Endo, Fujisawa, 252-0882, Japan, department of Pathology and Biology of 
Diseases, Graduate School of Medicine, Kyoto University, Sakyo-ku, Kyoto 606-8501 , Japan. 

Cytoplasm contains a large number of macromolecules at extremely high densities. How this striking nature 
of intracellular milieu called macromolecular crowding affects intracellular signaling remains 
uncharacterized. Here, we examined the effect of macromolecular crowding on ERK MAPK 
phosphorylation by MEK MAPKK. Addition of polyethylene glycol-6000 (PEG-6000) as a crowder to mimic 
intracellular environments, elicited a biphasic response to the overall ERK phosphorylation rate. 
Furthermore, probability of processive phosphorylation (processivity) of tyrosine and threonine residues 
within the activation loop on ERK increased non-linearly for increasing PEG-6000 concentration. Based on 
the experimental data, we developed for the first time a mathematical model integrating all of the effects of 
thermodynamic activity, viscosity, and processivity in crowded media, and found that ERK phosphorylation 
is transition-state-limited reaction. The mathematical model allows accurate estimation of the effects of 
macromolecular crowding on a wide range of reaction kinetics, from transition-state-limited to 
diffusion-limited reactions. 



Cells sense a variety of extracellular signals by receptors, and input signals to a system of intracellular 
reaction network, e. , an intracellular signal transduction cascade. These input signals are processed by the 
intracellular signal transduction cascades to drive the cells to exhibit specific phenotypes. A collapse of the 
system by gene mutations leads to pathological outcomes, including autoimmune disease and cancer 1 . Therefore, 
a better understanding of the intracellular signal transduction cascade is one of the critical issues for the control of 
diseases. The intracellular signal transduction cascade consists of a chain of diffusion, protein-protein interaction 
and enzymatic reactions, and thus it can be described by ordinary and partial differential equations, and solved 
numerically in order to analyze its dynamics 2 4 . Such a systems biology approach with experimentally- verified 
kinetic parameters provides a basis for understanding the signal transduction cascade 5 8 . 

The cytoplasm is highly crowded with macromolecules, such as proteins, lipids, nucleic acids and so on 9,10 . The 
concentration of proteins in the cytoplasm is over 100 mg/ml, with macromolecules making up 30-40% of the 
cytoplasmic volume 11,12 . This congested condition, often referred to as macromolecular crowding, results in 
enzymatic reaction kinetics that differ significantly from those described previously in ideal-dilute solution. 
For example, macromolecular crowding affects the enzymatic reaction rate through two factors: thermodynamic 
activity and viscosity 1316 . The former factor, thermodynamic activity, compensates for the concentration in an 
ideal solution to yield an effective concentration in a real solution 1517 . Macromolecular crowding increases 
thermodynamic activity mainly through the excluded volume effect, and consequently leads to an increase in 
a reaction rate. The latter factor, viscosity, reduces an encounter rate between the enzyme and the substrate 15 ' 18 . 
Therefore, the reaction rate is negatively correlated with viscosity. 

Processive reactions can be observed in phosphorylation 19-20 , polyubiquitination 21 , and so on. These multisite 
posttranslational modifications of signaling proteins are a common mechanism by which a complex cellular 
phenotype can be exhibited 22 . We demonstrated by a computer simulation that crowding environment increased 
the processivity in two-step phosphorylation reactions of ERK MAP kinase 23 . This prediction was validated by 
our experimental results 24 . 

Although macromolecular crowding has been suggested to have a great influence on the intracellular signaling 
pathway, it is not fully understood how the macromolecular crowding affects the overall reaction rate and 
processivity of enzymes. Here, we propose mathematical models of reaction kinetics that account for the effect 
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of thermodynamic activity, viscosity and processivity in an envir- 
onment with macromolecular crowding, and provide empirical vali- 
dations by employing in vitro ERK MAP kinase phosphorylation. 

Results 

Quasi-processive phosphorylation of ERK under macromolecular 
crowding. We have recently shown that Tyr and Thr residues within 
the activation loop of ERK MAP kinase were phosphorylated by 
MEK in a quasi-processive manner in mammalian cells 24 . In the 
first step, MEK binds to and phosphorylates ERK to generate tyro- 
sine monophosphorylated ERK (pY-ERK), and then dissociates from 
pY-ERK (Fig. 1A). Under macromolecular crowding, diffusion of 
MEK and pY-ERK is restricted. Therefore, rebinding of the MEK 
and the pY-ERK, which will not substantially happen in ideal solutes, 
occurs with some probability, and consequently the MEK proces- 
sively phosphorylates pY-ERK, producing the tyrosine and 
threonine bisphosphorylated ERK (pTpY-ERK). 

We analyzed the effects of macromolecular crowding by simplified 
quasi-processive phosphorylation model as follows. First, the rate 
constants of two successive phosphorylation reactions of ERK in 
the absence of a crowder agent were defined as k 1 and k 2 , respectively 
(Fig. IB, left). Our previous analysis demonstrated that these rate 
constants were first-order parameters, since the Michaelis constants 
in these reactions were significantly higher than the concentration of 
substrate, np-ERK, used in this study 24 . Under this condition, the kj 
and k 2 values can be determined by fitting experimental data (Fig. IB, 
left). Second, we built a simplified quasi-processive model, taking the 
crowding factor (c) and processivity (p) into consideration; c was 
defined as a dimension-less value that compensates for the increase 
in or decrease of reaction rates under macromolecular crowding in 
comparison with the reaction rates in the ideal-dilute solution, c 
reflects the effects of both thermodynamic activity and viscosity on 
the reaction rate 18 . Therefore, the rate constants in all reactions 
should be multiplied by c under a crowding condition (Fig. IB, right). 




c : Crowding factor 
ki, kz. phosphorylation rate p : Processivity (0 < p < 1 ) 



Figure 1 | Quasi-processive phosphorylation model under the condition 
of macromolecular crowding. (A) MEK (white) binds to and 
phosphorylates np-ERK (blue) at the tyrosine residue in the activation 
loop, and then dissociates from the product, pY-ERK (green). Under the 
condition of macromolecular crowding, the diffusions of both MEK and 
pY-ERK are restricted by the cage effect of the crowder molecules (grey). 
Thus, MEK rebinds to and phosphorylates pY-ERK with high probability. 
Consequently, MEK phosphorylates ERK in a processive manner, and 
produces pTpY-ERK (red). (B) The distributive (left) and quasi-processive 
phosphorylation model (right) are represented with parameters, including 
the phosphorylation rate constants, fcj and k 2 , crowding factor, c, and 
processivity, p. 



p is defined as the probability of a processive reaction; i.e., the pro- 
bability that pTpY-ERK will be generated directly from np-ERK 
(Fig. IB, right). The second reaction of processive phosphorylation 
would be on the order of 10-0.1 msec 23 , and much faster than kj and 
k 2 . For this reason, the reaction rate of the processive reaction can be 
described as k*p. According to this definition, the rate constant of 
pY-ERK generation must be multiplied by (1-p). With the k 1 and k 2 
values determined in the absence of a crowder, c and p values can be 
obtained by fitting the experimental data with the quasi-processive 
model 24 . 

Measurements of crowding factor and processivity by in vitro 
ERK phosphorylation in the absence or presence of a crowder, 
poly-ethylene glycol 6000. To evaluate the effect of macromole- 
cular crowding on c and p, the kinetics of ERK phosphorylation in 
vitro were studied in the absence or presence of a crowder agent, 
poly-ethylene glycol 6000 (PEG-6000) 24 . A constitutively active 
mutant of MEK1 (MEK) and a kinase-dead mutant of GST-fused 
ERK2 (GST-ERK) were used as the kinase and substrate of the in 
vitro phosphorylation reaction, respectively. Samples were subjected 
to phospho-affinity gel electrophoresis with Phos-tag compound, 
which clearly separates four phospho-isoforms of ERK (Fig. 2A). 
The fractions of phospho-isoforms of ERK, namely, non-phosphorylated 
ERK (np-ERK), pY-ERK and pTpY-ERK, were quantified by an 
Odyssey Infra-red Imaging System (Fig. 2B, and Supplementary 
Fig. SI online). Under our experimental condition, the order of 
electrophoretic mobility in the Phos-tag gel was, from the slowest 
to fastest, pY-ERK, pTpY-ERK, pT-ERK (if any), and np-ERK 25 - 26 . 
Two phosphorylation rates, k : and k 2 , were determined from the 
data obtained in the absence of PEG-6000 (Fig. 2B, left), c and p at 
each concentration of PEG-6000 were determined by fitting the ex- 
perimental data with the simplified quasi-processive model (Fig. IB, 
right and Supplementary Fig. SI and Dataset online). 

PEG-6000 showed a dose-dependent biphasic effect on c with the 
highest acceleration at approximately 10% (v/v) (Fig. 2C, green dots). 
This result is in agreement with the scheme documented by Minton 
et al. 18 . p increased steeply with increasing PEG concentration 
(Fig. 2C, red dots). 

Formulation of the crowding effect on the overall reaction rate. 

We first formulated the effect of crowding on the overall reaction 
rates, namely, c in this study. Minton and coworkers have established 
a basis for the mathematical framework of macromolecular 
crowding. Here, a biphasic effect on the enzymatic reaction rate is 
explained by increasing excluded volume and viscosity 18 . To 
formulate them, a simple enzymatic reaction is considered: 

E + S^ES — >£ + P, 

where E, S, P and ES represent the enzyme, substrate, product and 
enzyme-substrate transition state complex. 

When the encounter rate between E and S is much higher than the 
conversion rate of ES to E + P, the reaction is said to be transition- 
state-limited. Under such a condition, the overall rate constant, kf, is 
obtained by 

d -B= kf =k ts =klr, (i) 

where k ts is the transition-state-limited rate constant, and k° ts is the 
limiting value of k ts in the absence of an added crowder 13,27 . T is 
defined by 



where y B , y s , and y ES are the activity coefficients of the enzyme, sub- 
strate and enzyme-substrate transition state complex, respectively 27 . 
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Figure 2 | In vitro ERK phosphorylation by MEK in the presence of PEG-6000 used as a crowder agent. (A) ERK was phosphorylated by MEK at the 
indicated concentration of PEG-6000 (v/v). The samples were subjected to Western blotting analysis with Phos-tag. (B) The fractions of phospho- 
isoforms of ERK determined in (A) are plotted against time. Lines represent the fractions of phospho-isoforms of ERK obtained by fitting the 
experimental data with the quasi-processive model in Fig. IB. Raw and fitted data correspond to Exp #2 in Supplementary Dataset. (C) The crowding 
factor (c, green) and processivity (p, red) at each PEG-6000 concentration were determined by fitting as described in the text. Additional data are included 
in Supplementary Figure SI. Dashed lines are drawn for emphasis. 



The crowder and enzyme are considered as hard spherical particles 
with radius r c , and r E . The activity coefficient of the enzyme in a fluid 
of the crowder occupying fraction (p of total fluid volume is given by 



where 



\ny E = - [n(l-<p)+A l Q+A 2 Q 1 +A } Q\ 



Q=?/{l-9), 
A 1 =R 3 + 4.5R 2 + 6R+1.5, 
A 2 = 3R 3 +9R 2 + 4.5R, 
A 3 = 3R 3 +4.5R 2 , 



and 



R = r E /r C - 



(3) 



(4) 



(5) 



Similarly, the activity coefficients of the substrate and enzyme-sub- 
strate complex, y s and y ES , are also obtained with their radii, r s and 
r ES . We calculated the T value in the presence of crowder agents of 
different sizes, 5 kD, 20 kD, and 100 kD (Fig. 3A). For this calcula- 
tion, we assumed that all molecules were rigid spherical particles, in 
which radii were proportional to the cube root of their molecular 
weights. Consistent with previous reports 27 , the F value increased 
monotonously with increasing volume fraction of crowder agent and 
correlated inversely with the molecular weight of the crowder 17 . 

When the encounter rate between the enzyme and substrate is 
small relative to the conversion rate of enzyme-substrate complex 
to product, the reaction is said to be diffusion-limited. Under this 
condition, the encounter rate, k enc , which is proportional to the sum 
of the diffusion coefficients of enzyme and substrate, dictates the rate 
of reaction. The diffusion coefficient is inversely proportional to the 
viscosity of the solvent according to the Stokes-Einstein relationship. 
As an empirical approximation 28,29 , we may describe the rate of reac- 
tion with an exponential decay of crowder concentration 

kf = k em ~k° enc exp(-g<p), (6) 



where k° enc is the encounter rate constant in the absence of a volume- 
excluding background species, and g is a constant value that is a 
function of the relative sizes and shapes of the enzyme, substrate 
and crowder molecule 27 . We quantified the relative viscosity of 
PEG-6000 to water by an Ostwald viscometer, and obtained a £ value 
of 0.11 by fitting with the exponential function (Fig. 3B). Based on 
Eq. 6, the overall rate constant decreased exponentially with an in- 
creasing fraction of PEG volume under a diffusion-limited condition. 

Eqs. 1 and 6 did not recapitulate the biphasic effect of the overall 
reaction rate independently of each other. To the best of our know- 
ledge, the biphasic effects as observed in the experimental results 
have not been explicitly described by a mathematical framework, 
though such opposing effects of thermodynamic activity and viscos- 
ity on reaction rate have been well documented 1518,28 . We considered 
an intermediate condition between a transition-state-limited reac- 
tion and a diffusion-limited reaction. Noyes 30 indicates the relation- 
ship between the conversion and encounter rates as follows: 



K f — K ts + K enc ' 



k f = 



h Ic 

n -ts n -enc 

kts + kern 



Combining Eqs. 1, 6 and 7 gives the following equation: 



(7) 



(8) 



Using the overall reaction rate k°j and ratio of the encounter rate to 
conversion rate 0 , Eq. 8 is written as 



(l + 0)rexp(- 



fc° 



where 



kf r + 0exp(-g(p) 



^0 _ Venc 



(9) 



(10) 
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Figure 3 | Effects of macromolecular crowding on the overall reaction 
rate. (A) The T value that represents the thermodynamic activity was 
numerically simulated, and semi-log plotted against the crowder 
concentration according to Eq. 2-5. The parameters of the numerical 
simulations were as follows: r MEK = 45 1 ' 3 , r E RK — 42 1 ' 3 , Tmek-erk = 87 1 ' 3 , 
reader = 5 1/3 (blue), 20 1 ' 3 (green), 100 1 ' 3 (red). (B) The relative viscosity of 
PEG-6000 to water was measured by an Ostwald viscometer, and semi-log 
plotted against the PEG-6000 concentration as dots with s.d. (N = 3). The 
red line is fitted to the experimental data with a single exponential function 
(slope = 0.11). (C) The crowding factor, c, was calculated according to Eq. 
12 as a function of PEG-6000 concentration. Red and blue lines represent 
the effects of thermodynamic activity (8^ 0 in Eq. 12) and viscosity 
(0<g; 0) on the reaction rate, respectively. The green bold line is fitted to the 
experimental data (green dots as in Fig. 2C) with a 8 value of 2.26. Relative 
reaction rates calculated at the indicated 8 value are plotted as green dashed 
lines. The parameters of the numerical simulations are as follows: 
r MEK = 45 1 ' 3 , r ERK = 42 1 ' 3 , tmek-erk = 87 1 ' 3 , r croder = 6" 3 , g = 0.11. 

and 

(11) 

Based on the definition of c, we derived an equation of c as follows: 

_k f _ {\ + o)TtM-g<p) ri2l 

kj r+0ex P (-g<p) { ' 

Importantly, under the transition-state-limited reaction (0s> 1), Eq. 
9 corresponds to Eq. 1, whereas under the diffusion-limited reaction 



(#<C 1) Eq. 9 gives Eq. 6. According to the value of T and the g value 
in PEG-6000 solution, we calculated the c value at different concen- 
trations of the crowder with a changing 0 value as a variable (Fig. 3C, 
green lines). These results provided reasonable fits to the effect of 
crowder on the overall reaction rate in both the transition-state- 
limited and diffusion-limited reaction. By fitting the experimental 
data shown in Fig. 2C to this model, we obtained a 0 value of 2.26 
(Fig. 3C, bold green line). The value implicated the transition-state- 
limited reaction to some extent in ERK phosphorylation by MEK. 

Formulation of processivity. We have defined the processivity as a 
probability of two-step successive phosphorylation reactions from 
np-ERK to pY-ERK to pTpY-ERK by an identical MEK. To 
formulate the processivity, we consider the situation where ERK 
has just been phosphorylated by MEK and converted to pY-ERK. 
Particles of MEK and pY-ERK are treated as point-like particles. 
These particles diffuse normally at all times in this model. Since 
MEK contains only a single ATP-binding pocket 31 , after the 
phosphorylation reaction, MEK must release ADP and bind to 
ATP in a manner dependent on an ordered bi-bi reaction. We 
define the period required for the replacement of ADP with ATP 
as a relaxation time, T reJ . The three-dimensional probability 
distribution of pY-ERK after T rei is then given by 

p ™ ( -- H (^^ ex K-4^)' (i3) 

where r is the distance from MEK and D is the sum of the diffusion 
coefficient of MEK and pY-ERK in crowder solvent (Fig. 4A). We 
assume a reaction radius, r reac , in which MEK rebinds to and 



A 




Relaxation time 

t = Ira/ 




0 5 10 15 20 25 
PEG-6000 concentration (%) 



Figure 4 | Effect of macromolecular crowding on processivity. 

(A) Schematic view of the processive phosphorylation. (Left) MEK 
phosphorylates the tyrosine residue of ERK at t = 0. (Right) The relaxation 
time, T re z, is defined as the period during which MEK releases ADP and 
binds to ATP. pY-ERK dissociated from MEK moves away by the diffusion 
coefficient of D. At t = T re /, MEK processively phosphorylates pY-ERK, if 
pY-ERK is located within the distance of r rean the reaction radius. 

(B) According to Eq. 20, processivity was calculated as a function of 
PEG-6000 concentration. The red bold line is fitted to the experimental 
date (red dots as in Fig. 2C) with a T re ; value of 0.020 msec. Processivity 
values calculated at the indicated x re i values are depicted as red dashed lines. 
The parameters of the numerical simulations are as follows: r rmc =10 nm, 
D wMer = 10 um 2 /sec, g = 0.11. 
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processively phosphorylates pY-ERK within a time of T re; . The 
processivity, p, is then given by 

2n x r reac 

pYERK(D,r,z re ,)-r 2 sin OdfydOdr, (14) 



p{P,r reac ,z rel ) = 

piDSnacTrei) = Erf 



0 0 0 
' rrur 



\2^Dz re [) yj TlDz re l 

where Erf is an error function given by 



exp 



f r 2 \ 

I reac \ 



Erf(x) 



exp( — s 2 \ds. 



(15) 



(16) 



Taylor series of an error function and exponential function are given 
by 

2 / x 3 x 5 x 7 



Erf(x) = —=\x 1 h 

Jy ' V 3 10 42 



exp(x) = V - 

z — j n 

n = 0 



00 „n 



17 



(18) 



Based on Eqs. 17 and 18 for small x< 1, the Eq. 15 becomes 



?(AWfti)' 



4,/tiDV. 



^T,D-^. (19) 



The relations in Eq. 19 were consistent with a part of the results 
reported previously 23 ; the rebinding dynamics of MEK and pY- 
ERK is derived from a three-dimensional random walker returning 
to the origins. Because we simplified MEK and pY-ERK as point 
sources, we could not obtain the V in decay dynamics of processivity 
in a range of times shorter than the time to travel a molecular dia- 
meter 23 . The diffusion coefficient is inversely proportional to the 
viscosity of the solvent according to the Stokes-Einstein relationship. 
Thus, Eq. 15 becomes 



p(.D water ,r rmc ,T rd ,(p) =Erf 



Dwater^rel J 



exp(gffl) 



(20) 



exp 



r Lc e Mg<p) \ 

^E) wa t er T re l J 



where D water is the sum of the diffusion coefficient of MEK and pY- 
ERK in water. To evaluate this model, we calculated processivity 
according to Eq. 20. D water was set at 10 um 2 /sec based on the vis- 
cosity of water and the size of MEK and ERK. In this model, we did 
not consider the effect of viscosity on diffusion coefficient of ADP/ 
ATP, i.e. z re i, because of their high concentration and fast diffusion 
(see Discussion). We set the reaction radius r reac at 10 nm, which was 
comparable to the sum of the radius of MEK and ERK molecules. 
Processivity was obtained by changing the value of z reh namely, the 
reactivation time of MEK, and plotted against the concentration of 
PEG-6000, which increased viscosity (Fig. 4B). The processivity ele- 
vated with increasing concentration of PEG-6000. Furthermore, the 
model demonstrated a best fit to the experimental data at the z re i 
value of 0.020 msec (Fig. 4B, bold line). 

A full reaction kinetics model of ERK phosphorylation. In the 

quasi-processive phosphorylation model described in this study, 
the association and dissociation of ERK to MEK were neglected for 
the simplicity. Therefore, we performed similar analyses based on the 
full reaction kinetics of ERK phosphorylation model. Kinetic para- 
meters were fitted with experimental data in the absence of PEG- 
6000 (Fig. S2A, left), followed by the fitting of crowding factor, c, and 
dissociation constant, k b -, in the full model with experimental data 



under the condition of PEG-6000 (Fig. S2A right). Of note, the 
reaction model can strictly define processivity as k 2 l(k b - + k 2 ) with- 
out any need to introduce an empirical parameter of processivity in 
the simplified model (Fig. IB). The full model recapitulated time 
course of ERK phosphorylation in the absence or presence of PEG- 
6000 (Fig. S2B). The fitted values of crowding factor (Fig. S2C) and 
processivity (Fig. S2D) provided consistent values of 0 (1.80) and 
z rd (0.024 msec) in Eq. 12 and Eq. 20, respectively. 

Analysis of crowding factor and processivity based on the scaled 
particle theory. The crowding effect (Eq. 12) and processivity (Eq. 
20) explicitly include viscosity, which is macroscopic parameter 
described by the Stokes-Einstein law. Therefore, we next examined 
the microscopic effect of crowder on crowding factor and proces- 
sivity by estimating diffusion coefficients based on the scaled particle 
theory (SPT) 32 . Let r f represent the radius of the tracer species, Ar' the 
approximated length of the Brownian displacement 32 . The scaled 
particle theory yields 

ln(D/D water ) = -Ar' 

[A i +A 2 (2r t + Ar') + A 3 (3r 2 + 3r ( Ar' + Ar' 2 ) ] , 

where D and D water are the diffusion coefficient of tracer species in 
the presence or absence of background crowder. A, is given by 



Ao = -ln(l-5,), 
A 1 =6S 2 /(1-S 3 ), 
A 2 = [12Si/(l -%)] + [18S 2 /(1 -S3) 2 ] , 
A 3 = [8S 0 /(1 -S 3 )] + [24S!S 2 /(1 -S3) 2 ] + [2AS\/{\ -S 3 ) 3 ] , 

_ CbackNa 

V ~ 1000M' 



(22) 
(23) 
(24) 
(25) 
(26) 

(27) 



where r b is radius of background molecule, M is molecular weight, N a 
is Avogadro's number, c back is a concentration of background mole- 
cules [g/litter], and v is in cm -3 . We obtained almost same result of 
BSA self- diffusion as previously reported in Muramatsu and Minton, 
indicating the valid calculation (Fig. S3 A). Next, logarithmic ratio of 
D to D water of relative diffusion coefficient of MEK and ERK with the 
increase in PEG-6000 concentration was estimated by SPT (Fig. 
S3B). The microscopic model with SPT significantly underesti- 
mated the relative diffusion coefficient in the presence of PEG- 
6000 (Fig. S3B, red line) in comparison to the relative diffusion 
coefficients calculated by the experimental data with Ostwald visco- 
meter (Fig. S3B, blue line). We recalculated and fitted the crowding 
factor (Fig. S3C) and processivity (Fig. S3D) using SPT with the ex- 
perimental data; however there existed significant discrepancy be- 
tween the experimental data and mathematical model with the scaled 
particle theory. This could be due to steeper reduction of the relative 
diffusion coefficient between MEK and ERK obtained by SPT. From 
these results, we concluded, at this time, that the increase in viscosity 
under the condition of macromolecular crowding was the respon- 
sible factor for the crowding factor and processivity. The clear dissec- 
tion of difference between microscopic and macroscopic effect of 
molecular crowding is one of challenging tasks for future studies. 
Simulation of the enzymatic reaction at the particle level will enable 
to directly address this issue 23 . 

Analytical solution of the simplified quasi-processive phosphory- 
lation model. Finally, based on Eq. 12 and 20, analytical solutions of 
np-ERK, pY-ERK and pTpY-ERK concentration, [npERK], [pYE- 
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RK], and [pTpYERK], respectively, in the quasi-processive model 
(Fig. IB) are given by 

[npERK]{t)=ERK a -e~ cklt (28) 

e -c(k l+ h)t f e c-k 2 t _ e c-ht\ k r p _ l \ 
\p YERK] ( t) = ERK 0 ^— (29) 

\pTpYERK] (r) = ERK 0 

e -c(k, +k 2 )t^(h +fe)t( fcl + kl) + e ok,t ki (p _ i) + gc-fe* (jt 2 _ } (30) 

where the [npERK], [pYERK], and [pTpYEMC] at time = 0 are ERK 0 , 
0, and 0, respectively. k t and k 2 were first-order successive 
phosphorylation rates [/sec] of ERK by MEK based on our 
previous study 24 . To simplify these equations, we omitted the 
factor of MEK concentration, and included it implicitly in k 1 and 
k 2 . In general, analytical solutions of two successive reactions were 
obtained at c = 1 and p = 0 as follows: 



[npERK](t)=ERK 0 -e 



-hi 



(31) 



e -(h+k 2 )t/_ e k 2 t _|_ g fcit\ h 
\p YERK] (t) = ERKq — £ j (32) 

fci — K.2 

\pTpYERK}{t)=ERK 0 

e -(ti+fa)*{e(fe+fe)f(jfc 1+ fe 2 )_e*ifjfc 1 + fc 2e fe*}(33) 
fci-fc 2 

In the quasi-processive phosphorylation model, we assumed that the 
reaction rate of processive phosphorylation is equivalent to fcj*p, 
because the second reaction of processive phosphorylation would 
be much faster than fc ; and k 2 . These equations provided a good 
agreement with the dynamics of ERK phosphorylation at different 



PEG-6000 0% 



PEG-6000 6% 




10 20 
Time (min) 

PEG-6000 12% 



10 20 
Time (min) 

PEG-6000 18% 




np-ERK 
pY-ERK 
pTpY-ERK 



10 20 
Time (min) 



10 20 
Time (min) 



Figure 5 | Analytical solution of ERK phosphorylation by MEK under the 
condition of macromolecular crowding. The fractions of np-ERK (blue), 
pY-ERK (green) and pTpY-ERK (red) were calculated according to 
Eq. 21-23. Experimental data are the same in Figure 2A and 2B (Exp. #2 in 
Supplementary Data set). The parameters of the numerical simulations 
were as follows: r MEK = 45 1 ' 3 , r EliK = 42 1 ' 3 , r MEK _ ERK = 87" 3 , r croder = 6 m , 
g = 0.11, 6 = 2.26, r reac = 10 nm, D water = 10 unrVsec, T re j = 0.020 msec. 



concentrations of PEG-6000 (Fig. 5), supporting the validity of the 
model. 



Discussion 

In this study, we provided for the first time the mathematical frame- 
work for the rate constant and processivity at intermediate condi- 
tions between transition-state-limited and diffusion-limited 
reactions under the condition of molecular crowding. Moreover, 
we derived an analytical expression for the kinetics in ERK processive 
phosphorylation, demonstrating time courses of phosphorylation as 
a function of crowder concentration, <j>. These models recapitulated 
reaction kinetics of dual phosphorylation of ERK by MEK with or 
without crowder agent. However, slight, but not ignorable, disagree- 
ment is also evident between the experimental data and the model, 
especially in late phase of higher concentrations of PEG-6000, 12% 
and 18%. This result suggests that the PEG-6000 solution demon- 
strates a spatially heterogeneous environment, which entails the 
introduction of a time-dependent rate coefficient with a fractal kin- 
etics model 13,33 . In a previous simulation study, diffusion obstacles 
gave rise to anomalous diffusion 34 . Consistent with this idea, anom- 
alous diffusion in the cytoplasm was also observed by fluorescence 
correlation spectroscopy 35 37 . Moreover, it has recently been indi- 
cated that fractal structures exist in the nucleus 38,39 . These findings 
should pave the way for future work that will improve the model 
more quantitatively. 

We demonstrated a clear dependency of processivity in Eq. 19 on 
diffusion coefficient, D, reactivation time, T re ;, and reactivation 
radius, r reac . Membrane-anchoring or transmembrane proteins are 
well-known to show slower diffusion coefficients in comparison to 
those of cytoplasmic proteins 40 . Of note, many of receptor tyrosine 
kinases (RTKs) such as epidermal growth factor receptor form 
homo- and hetero-dimer upon ligand stimulation on membrane, 
and transphosphorylates multiple sites of tyrosine residues in 
RTKs 41 , strongly suggesting the possibility of processive phosphor- 
ylation in the RTK activation. Diffusion-controlled ligand-receptor 
kinetics supports this possibility 42 . Furthermore, rebinding event 
would take place in protein translocation on nucleic acid 43,44 . 
Intriguingly, by analogy to the diffusion coefficient, clustering or 
confinement of molecules in small compartment also induced the 
increase in processivity 45 . 

Taking the high concentration and fast diffusion rate of intracel- 
lular ATP molecules into consideration, the rate-limiting step of 
ADP/ATP exchange, i.e. reactivation time, T re ;, is likely to be deter- 
mined by the ADP release from MEK. The dissociation rate con- 
stants of ADP from kinases have been reported to vary over a 
wide range with an order of magnitude of 3 to 4 (0.23 sec -1 to 
140 sec- 1 ) 46 - 50 . The ADP/ATP exchange rate predicted in this study, 
50000 sec -1 , which is the reciprocal of z re i is much larger than the 
dissociation rate constants reported previously. There are at least two 
possible explanations for this discrepancy. First, the reaction radius, 
r reaa ma y have been underestimated in this study. Currently, we do 
not have any methods to quantitatively estimate the reaction radius. 
Second, the phosphorylation reaction model may not be correct. We 
employed an ordered bi-bi sequential reaction of processive phos- 
phorylation, in which MEK dissociated from pY-ERK, followed by 
exchange of ADP with ATP to phosphorylate pY-ERK. If MEK pro- 
cessively phosphorylates ERK through a random bi-bi mechanism, the 
processivity will increase even under a slow ADP/ATP exchange rate. 
In fact, Horiuchi et al. have indicated that the mechanism of MEK 
turnover involves random addition of substrate and ATP 51 . Although 
they have also suggested the ordered release of product from MEK 51 , 
effects of crowding on the mechanism of MEK turnover have not been 
fully characterized, and should be considered in a future study. 

The mathematical models developed here integrate three effects of 
molecular crowding on the reaction rate — namely, thermodynamic 
activity, viscosity and processivity. These models were consistent 
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with experimental data supporting the processive ERK phosphoryla- 
tion by MEK. Our findings will contribute to the effort to quantita- 
tively simulate intracellular signal transduction. 

Methods 

Plasmids. pGEX4T3-ERK2KR, a kinase -deficient mutant (K57R) of Xenopus ERK2, 
and pET- 12HisFLAG-MEKlSDSE, a constitutively active mutant (S218D/S222E) of 
Xenopus MEK1, have been purified as described previously 24 . 

Reagents and antibodies. Phos-tag acrylamide was obtained from the Phos-tag 
Consortium (http://www.phos-tag.com/). Fifty percent polyethylene glycol 6000 
(PEG-6000) solution was purchased from Hampton Research (Aliso Viejo, CA). 
Anti-ERKl/2 and anti-phospho-ERKl/2 {Thr202/Tyr204 and Thrl85/Tyrl87, 
respectively) were purchased from Cell Signaling Technology (Beverly, MA). 
LI-COR-blocking buffer and IRDye680- and IRDye8 00 -conjugated anti-rabbit and 
anti-mouse IgG secondary antibodies were obtained from LI-COR Bioscience 
(Lincoln, NE). 

Purification of recombinant proteins and in vitro kinase assay. GST-ERK2KR and 
12His-MEKlSDSE were prepared as described previously 24 . For the in vitro kinase 
assay, 1 uM GST-ERK2, 0.1 uM 12HisFLAG-MEKlSDSE and different 
concentrations of PEG-6000 (0—20% v/v) were incubated in in vitro kinase buffer 
(50 mM Tris [pH 7.5], 10 inM MgCl 2l 0.02% BSA 0.2 mM DTT). The reaction was 
started by adding 5 X ATP solution (5 mM ATP, 50 mM MgCl 2 , 25 mM Tris 
[pH 7.2], 0.15 M NaCl) and stopped at the indicated time point by adding SDS 
sample buffer. 

Phos-tag polyacrylamide gel electrophoresis. Phos-tag polyacrylamide gel 
electrophoresis was performed essentially as described previously 24,25 . Briefly, 50 uM 
Phos-tag and 100 uM MnCl 2 were added to conventional SDS -polyacrylamide 
separation gels according to the manufacturer's protocol, respectively. Proteins were 
detected by using an Odyssey Infrared Imaging System (LI-COR). 

Measurements of viscosity. Viscosities of the aqueous PEG-6000 solutions relative to 
water were determined by using an Ostwald viscometer at 37°C in a water bath. 

Parameter search. Numerical simulation and parameter search were implemented by 
MATLAB software with odel5s function and fminsearch function, respectively. 
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